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Abstract 

Cells  were  life-cycled  cells  using  profiles  with  a  3,  6,  or  9%  change  in  state  of  charge  (ASOC)  at  40,  50,  60,  and  70  °C.  From  the  voltage 
response  of  the  cells  to  the  life-cycle  profile  at  each  temperature,  we  separated  the  overall  impedance  rise  into  two  simpler  terms,  R0  (ohmic) 
and  Rp  (polarization),  using  an  equivalent  circuit  model.  The  R0  data  tend  to  follow  the  expected  trends  (40  >  50  >  60  >  70  °C).  Although 
the  Rp  data  trends  show  that  Rp  can  either  decrease  or  increase  asymptotically  with  time,  the  overall  temperature-dependent  behavior  is  similar 
to  that  of  R0.  We  illustrate  the  types  of  processes  that  can  occur  in  one  lithium-ion  cell  chemistry.  Based  on  the  initial  rates,  the  processes  are 
complex.  The  R0  term  dominates  the  observable  cell  impedance,  but  Rp  adds  a  non-trivial  contribution. 
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1.  Introduction 

Since  their  introduction  in  1991  [1],  Li-ion  batteries  have 
found  uses  in  many  consumer  electronic  applications.  The 
batteries  generally  consist  of  a  metal  oxide  cathode,  a  liquid 
organic  solvent  containing  a  dissolved  ionic  salt,  and  a 
carbon-based  anode.  The  US  Department  of  Energy — sup¬ 
ported  Partnership  for  the  New  Generation  of  Vehicles 
(PNGV)  and  others  have  been  evaluating  this  class  of 
batteries  for  use  in  hybrid-electric  vehicles.  For  this  applica¬ 
tion,  the  batteries  must  have  a  useful  life  of  15  years. 
However,  during  their  lives,  these  cells  exhibit  an  impedance 
increase.  Thus,  a  15-year  life  remains  a  challenge. 
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Three  national  laboratories,  Argonne  National  Laboratory 
(ANL),  Idaho  National  Engineering  and  Environmental 
Laboratory  (INEEL),  and  Sandia  National  Laboratories 
(SNL),  collaborated  to  evaluate  and  understand,  among 
other  phenomena,  the  causes  of  the  impedance  rise  in  these 
cells.  A  number  of  18650-sized  cells  containing  a  commer¬ 
cially  available  chemistry  were  subjected  to  accelerated 
cycle  and  calendar  life  testing  in  this  study.  The  papers 
describing  this  study  and  some  of  the  results  appear  else¬ 
where  [2,3].  These  studies  showed  that  the  discharge  impe- 
dance  increase  is  proportional  to  (time)  ,  regardless  of  the 
aging  test.  The  authors  infer  that  the  most  likely  source  of 
this  impedance  increase  is  the  growth  of  a  solid  electrolyte 
interface  (SEI)  layer  on  the  surface  of  one  or  both  of  the 
electrodes. 

The  PNGV  Battery  Test  Manual  [4]  describes  a  simple 
model  to  estimate  the  Ohmic  and  polarization  impedances 
that  constitute  the  observable  discharge  and  regeneration 
(regen)  resistances.  The  Lumped-Parameter  Model  (LPM) 
[5]  attempts  to  linearize  the  behavior  and  response  of  a 
battery  during  a  simple,  repeatable  current  profile  such  as 
that  used  in  the  hybrid  pulse  power  characterization  test  [4] . 
The  life-cycle  profiles  used  in  this  study  were  based  on 
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current,  not  power,  so  the  LPM  was  used  to  analyze  the  life- 
cycle  voltage  and  current  versus  time  data.  From  this 
analysis,  we  show  that  the  underlying  mechanisms  are 
sensitive  to  the  state  of  charge  (SOC),  amount  of  charge 
removed  during  the  cycling  (ASOC),  and  the  temperature. 

2.  Experimental 

Details  of  the  cell  chemistry  and  the  conditions  that  the 
cells  were  exposed  to  are  given  in  [2].  Briefly,  135  0.9  A  h, 
18650-sized  cells  were  built  to  ANL  specifications.  The 
cells  consisted  of  a  Li(Co0.2Ni0.8)O2  plus  SFG-6  +  carbon 
black  graphite  cathode,  1  M  LiPF6  in  EC:DEC  (1:1  by  wt) 
electrolyte,  a  Celgard  separator,  and  a  MCMB-6  +  SFG-6 
graphite  anode.  The  active  area  was  678  cm2.  The  cycle- 
life  cells  were  aged  at  60  and  80%  SOC  at  40,  50,  60,  and 
70  °C.  The  cells  were  cycled  using  three  profiles,  3,  6,  and 


9%  ASOC,  as  shown  in  Fig.  1.  Three  to  six  cells  were 
exposed  to  all  aging  conditions.  Every  4  weeks  (2  weeks  at 
70  °C),  the  cells  were  cooled  to  25  °C  to  conduct  reference 
performance  tests.  During  the  life-cycle  experiments,  data 
from  every  100th  profile  were  collected  at  high  rates  (ca.  10 
samples/s).  These  data  were  extracted  and  analyzed  in 
terms  of  the  LPM. 

The  equivalent  circuit  used  in  the  LPM  is  given  in  Fig.  2a, 
where  OCV  is  the  ideal  open  circuit  voltage;  R0  is  the 
internal  resistance  due  to  ohmic,  interfacial,  and  kinetic 
impedances  in  the  cell;  Rp  is  the  concentration  polarization 
of  the  electrolyte  and  electrodes;  C  is  a  capacitance  included 
to  simulate  the  time  constant  during  the  polarization  of  the 
electrolyte  and  electrodes;  /L  is  the  load  current;  Ip(t )  the 
current  through  the  polarization  resistance;  VL  load  voltage; 
Rl  load  resistance;  t  the  time;  and  1/OCV'  is  related  to  the 
changes  in  the  diffusion  coefficient  of  Li+  in  the  electrodes 
with  time. 
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Fig.  1.  Cycle-life  profiles  used:  (a)  3%  ASOC,  (b)  6%  ASOC,  (c)  9%  ASOC. 


154 


I.  Bloom  et  al. /  Journal  of  Power  Sources  111  (2002)  152-159 


1/OCV’ 


Fig.  2.  (a)  Equivalent  circuit  used  for  the  lumped-parameter  battery  model;  (b)  comparison  of  calculated  and  measured  cell  voltages. 


The  relationship  between  VL  and  most  of  the  parameters 
listed  is  described  by  Eq.  (1): 


yL  =  ocv  -  ocv' 


ILdt 


—  Rq  [/l]  —  ^p[^p 


altered  while  monitoring  r  ;  if  a  more  complex  form  of  the 
function  yielded  a  significantly  higher  r  (e.g.  0.85  versus 
0.95),  then  the  new  function  was  selected  as  the  one  more 
representative  of  the  process. 


The  parameters,  OCV,  OCV',  Rp,  R0 ,  and  x  were  calculated 
by  using  the  linear  regression  function  (LINEST)  in  Micro¬ 
soft  Excel.  An  iterative  computer  algorithm  of  successive 
approximations  (t  +  At)  was  written  to  estimate  x  (the 
initial  value  of  At  was  1).  In  this  method,  the  values  of 
successive  regression  coefficients,  r2,  were  monitored  as  a 
function  of  t.  As  long  as  r  increased,  the  value  of  At  was  not 
changed  and  was  added  to  t.  When  r  decreased,  the  value  of 
At  was  halved  and  negated.  This  process  was  continued  until 

O  _ 1  A 

two  successive  values  of  r  did  not  differ  by  more  than  10 
Typically,  the  final  values  of  r  were  0.95  or  better.  Fig.  2b 
shows  that  the  calculated  values  agree  well  with  the  mea¬ 
sured  cell  voltage  values. 

After  the  calculations  were  completed,  the  individual 
values  of  R0  and  Rp  were  averaged.  The  time-dependencies 
of  Rp  and  R0  were  determined  by  curve  fitting.  The  type  of 
equation  used  is  given  in  Eq.  (2): 

Q  =/(*),  (2) 

where  Q  =  R0  or  Rp,  and  f(t)  is  of  the  form  at  +  b , 
a\n(t)  +  b,  at1^2  +  b,  at1!1  +  b\n(t)  +  c  (< a ,  b,  c  =  con¬ 
stants,  t  =  cycle  count),  or  some  linear  combination  of 
them.  The  form  of  the  equation  was  selected  based  on  its 
value  of  r  resulting  from  fits  of  the  data.  The  form  of  Q  was 


3.  Results  and  discussions 

An  abbreviation  system  is  used  below  to  denote  the 
experimental  conditions.  It  signifies  the  SOC,  ASOC,  and 
temperature  of  the  experiment.  For  example,  “80360”  sig¬ 
nifies  80%  SOC,  3%  ASOC  at  60  °C,  and  60940  signifies 
60%  SOC,  9%  ASOC  at  40  °C.  If  the  temperature  field  is 
denoted  as  “XX,”  then  we  are  referring  to  the  entire  group  of 
cells  at  a  particular  SOC  and  ASOC.  It  should  be  noted  that 
data  from  the  first  ~5000  cycles  or  the  first  4- week  soak 
interval  (2  weeks  at  70  °C)  were  used  (approximately 
10,000-25,000  cycles).  After  that,  the  data  were  either 
obscured  by  experimental  artifacts  or  too  complex  to  ratio¬ 
nalize.  Hence,  the  analyses  indicate  the  initial  mechanisms 
and  rates  only.  For  clarity,  the  remaining  discussion  is 
divided  according  to  source  of  the  impedance,  RQ  or  Rp. 

3.1.  Ohmic  impedance  rise  (R0) 

Typical  plots  of  R0  versus  cycle  count  for  cells  cycled  at 
60%  SOC  are  given  in  Fig.  3;  plots  for  those  cycled  at  80% 
SOC,  in  Fig.  4.  Comparison  of  Figs.  3  and  4  shows  both 
similarities  and  differences.  The  general  shapes  of  R0  versus 
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Fig.  3.  Average  R0  vs.  cycle  count  for  cells  at  606XX  (XX  =  temperature  given  on  plot). 


cycle  count  are  similar,  but  their  dependence  on  temperature 
is  different.  In  Fig.  3,  the  60%  SOC  data  follow  a  simple, 
temperature-dependent  trend.  The  values  of  R0  decrease 
with  rising  temperature,  40  >  50  >  60  >  70  °C.  However, 
Fig.  4  shows  that  this  is  not  the  case  at  80%  SOC.  The  values 
of  R0  decrease  as  40  >  50  >  70  >  60  °C.  Exactly,  why  the 
80%  SOC  values  do  not  follow  the  expected  trend  is  not 
known. 

The  values  of  the  coefficients  of  the  fitting  equation, 
a\n(t)  +  bt1/2  +  c,  and  the  values  of  r2  from  curve  fitting, 
are  given  in  Table  1  for  cells  aged  at  60%  SOC.  The 
corresponding  data  from  the  80%  SOC  cells  using  the 
equation  at1!2  +  b  are  given  in  Table  2.  As  can  be  seen 
from  the  high  values  of  r  ,  the  fits  are  very  good.  A  plot 
showing  typical  the  agreement  between  the  data  and  the  fits 
is  given  in  Fig.  5. 

From  Tables  1  and  2,  we  can  see  that  the  form  of  the 
equation  that  describes  R0  differs  according  to  SOC;  the 
form  of  the  equation  for  60%  SOC  differs  from  that  for  80% 
SOC.  Though  the  values  of  the  coefficients  can  be  used  to 
compare  absolute  rates,  the  form  of  the  equation  reveals 
much  about  the  underlying  processes. 


Square-root-of-time-dependent  (parabolic  kinetics),  as 
well  as  logarithmic-with-time,  mechanisms  have  been  found 
for  the  growth  of  an  oxide  film  on  a  metal  surface  [6-14]. 
The  square-root-of-time-dependent  mechanism  tends  to 
have  a  strong  thermal  diffusion  component,  tends  to  follow 
Arrhenius-like  kinetics,  and  results  in  the  formation  of 
compact  scales.  In  the  logarithmic-with-time  case,  the  rate 
of  oxidation  starts  rapidly  and  soon  becomes  slow.  This  rate 
law  has  been  described  as  following  from  either  the  absorp¬ 
tion  of  reactive  species  or  the  effect  of  electric  fields  in  or 
across  the  oxide  layer. 

Applying  this  analogy  to  the  Li-ion  battery  case,  R0 
behaves  as  though  the  formation  of  a  thin  film  or  some 
other  resistive  layer  on  a  battery  electrode  may  occur  in 
much  the  same  manner  as  an  oxide  film  on  a  metal.  In  the 
60%  SOC  case,  the  form  of  the  equation  implies  that  the 
reaction  starts  very  quickly  and  then  slows  down.  The  initial 
reaction  between  the  electrolyte  and  electrode(s)  may 
involve  absorption  of  the  organic  electrolyte  and/or  LiPF6 
on  the  electrode  surface  as  the  rate-determining  step.  The 
absorbed  material  then  reacts  quickly  to  form  a  very  thin 
layer.  As  soon  as  sufficient  material  has  reacted — hence,  the 
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Fig.  4.  Average  R0  vs.  cycle  count  for  cells  at  803XX  (XX  =  temperature  given  on  plot). 
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Table  1 


Coefficients  for  R0  vs.  cycle  count  equations  for  60%  SOC 


SOC  (%) 

Temperature  (°C) 

ln(f) 

A2 

Constant 

r2 

3 

40 

5.84 

X 

i<r5 

5.86 

X 

icr5 

2.46 

X 

10“2 

0.97 

50 

3.77 

X 

10“4 

4.69 

X 

10“4 

2.07 

X 

10“2 

0.98 

60 

4.95 

X 

10“4 

2.24 

X 

icr5 

1.77 

X 

icr2 

0.91 

70 

-2.01 

X 

1(T4 

4.16 

X 

1(T5 

2.26 

X 

icr2 

0.97 

6 

40 

1.50 

X 

i(T3 

4.98 

X 

icr5 

2.41 

X 

icr2 

0.98 

50 

9.24 

X 

1(T4 

5.98 

X 

1(T5 

2.30 

X 

icr2 

0.99 

60 

3.63 

X 

1(T4 

6.16 

X 

1(T5 

2.42 

X 

icr2 

0.99 

70 

2.90 

X 

10“4 

4.18 

X 

10“5 

2.50 

X 

10“2 

0.99 

9 

40 

-6.74 

X 

1(T4 

1.80 

X 

1(T4 

3.24 

X 

icr2 

0.98 

50 

9.87 

X 

1<T5 

9.54 

X 

icr5 

2.83 

X 

icr2 

0.98 

60 

5.36 

X 

10“4 

5.93 

X 

10“5 

2.54 

X 

10“2 

0.96 

70 

3.81 

X 

10“5 

8.12 

X 

10“5 

2.52 

X 

10“2 

0.98 

Table  2 


Coefficients  for  R0  vs.  cycle  count  equations  for  80%  SOC 


SOC  (%) 

Temperature  (C) 

tm 

Constant 

2 

r 

3 

40 

1.11 

X 

10“4 

2.81  x  1(T2 

0.99 

50a 

6.20 

X 

10“5 

2.69  x  10“2 

0.91 

60 

5.69 

X 

icr5 

2.14  x  10“2 

0.96 

70 

9.07 

X 

1(T5 

2.46  x  1(T2 

0.92 

6 

40 

1.15 

X 

10“4 

3.52  x  1CT2 

0.96 

50 

1.49 

X 

1(T4 

3.17  x  1(T2 

0.92 

60 

9.11 

X 

1(T5 

2.64  x  1(T2 

1.00 

70 

1.36 

X 

10“4 

2.93  x  10“2 

0.91 

9 

40 

2.29 

X 

1(T4 

3.16  x  10“2 

0.91 

50 

1.27 

X 

icr4 

3.12  x  1(T2 

0.96 

60 

8.20 

X 

10“5 

2.95  x  10“2 

0.88 

70 

1.21 

X 

10“4 

2.64  x  10“2 

0.95 

aThis  condition  could  also  be  expressed  as  R0  =  1.70  x  10  4ln(t)  + 
1.70  x  10-2.  The  value  of  r2  is  0.91. 


layer  has  grown — the  reaction  slows  down.  At  this  point,  the 
diffusion  of  Li+  and  electrons  through  the  layer  become  the 
rate-limiting  step.  The  change  in  rate-limiting  step  is  reflected 
by  the  parabolic  term  in  the  equation.  Unfortunately,  our 
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Fig.  5.  Average  R0  vs.  cycle  count  and  fitting  curves  using  the  data  in 
Table  1  for  609XX  (XX  =  temperature  on  graph).  The  markers  represent 
the  values  of  the  average  R0  and  solid  lines  represent  the  fitting  curves. 


results  do  not  indicate  on  which  electrode  these  reactions 
occur.  Reports  in  the  literature  indicate  that  the  greatest 
impedance  rise  for  our  cell  chemistry  occurs  on  the  cathode 
due  to  SEI  layer  growth  [12]. 

From  the  80%  SOC  data,  the  form  of  the  equation 
indicates  that  parabolic  kinetics  are  operating.  That  is,  the 
diffusion  of  Li+  and/or  electrons  through  the  film  may  be  the 
rate-determining  step,  not  the  adsorption  of  materials. 
Equally  possible,  the  diffusion  of  electrolyte  solvent  and / 
or  salt  molecules  through  the  film  to  the  cathode  could  be  the 
rate-determining  step. 

Conceptually,  the  anode  could  also  play  an  important  role 
in  the  impedance  rise.  Similar  mechanisms  as  described  for 
the  cathode  may  be  operating.  Based  on  the  data,  we  cannot 
differentiate  which  electrode  contributes  the  most  to  the 
impedance  rise. 

Plotting  the  coefficients  of  the  t  ~  terms  on  Arrhenius 
plots  yields  additional  information  (Fig.  6).  At  60%  SOC, 
Arrhenius  kinetics  are  followed  at  3%  ASOC;  at  9% 
A  SOC,  Arrhenius  kinetics  are  followed  up  to  (and  including) 
60  °C.  Arrhenius  kinetics  are  not  followed  at  6%  ASOC 
(r2  =  0.11).  From  the  slopes  of  the  lines  in  Fig.  5,  the  thermal 


ASOC  at  60%  SOC. 
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activation  energies  (Ea)  for  R0  are  41.  3  ±  13.6  (r2  =  0.90) 
and  48.2  ±  3.2  kJ/mol  (r  =  1.00)  for  3  and  9%  ASOC, 
respectively.  These  values  are  the  same  within  experimental 
error. 

Note  that  the  coefficients  of  the  ln(7)  term  also  follow 
Arrhenius  kinetics.  From  the  slope  of  the  line,  the  value  of  Ea 
is  52.7  zb  8  kJ/mol  (r2  =  0.96). 

The  same  exercise  using  the  80%  SOC,  R0  data  shows 
that  only  9%  ASOC  seems  to  follow  Arrhenius  kinetics, 
and  then  only  up  to  60  °C.  Here,  the  value  of  Ea  for  R0  is 
44.6  zb  3.0  kJ/mol  (r2  =  1.00).  There  is  no  significant  dif¬ 
ference  between  this  and  the  previous  values  of  Ea. 

3.2.  Polarization  impedance  rise  (Rp) 

Typical  plots  of  Rp  versus  cycle  count  are  given  in  Figs.  7 
and  8  for  cells  cycled  at  60  and  80%  SOC,  respectively. 
Comparing  Figs.  7  and  8  shows  two  very  different  behaviors 
in  Rp.  In  Fig.  7,  the  value  of  Rp  decreases  asymptotically  with 
cycle  count;  in  Fig.  8,  there  is  a  slow  increase.  It  is  clear  that 
the  curves  in  Fig.  7  follow  an  expected  trend  based  on  tem¬ 
perature,  40  >  50  >  60  >  70  °C.  In  Fig.  8,  this  is  not  the  case. 

The  values  of  the  coefficient,  the  form  of  the  equation,  and 
the  values  of  r2  from  the  curve  fitting  are  given  in  Table  3  for 
cells  aged  at  60%  SOC.  As  explained  below,  there  are  two 
rows  of  coefficients  for  the  606XX  conditions,  one  for  an 
equilibration  process  and  the  other  for  a  degradation  process. 
The  corresponding  data  from  the  80%  SOC  cells  are  given  in 
Table  4.  As  can  be  seen  from  the  high  values  of  r  ,  the  fits  are 
very  good,  with  two  exceptions.  The  data  from  60340  seem 
constant  with  cycle  count,  and  the  data  at  80940  show  too 
much  scatter. 

From  the  data  in  Table  3,  the  60340  data  indicate  that  there 
is  no  dependence  of  Rp  on  cycle  count.  This  implies  that 
there  is  little,  if  any,  change  in  concentration  polarization  in 
the  cells.  As  the  temperature  increases,  the  amount  of 
concentration  polarization  increases,  most  likely  due  to 
higher  reaction  rates  for  the  decomposition  of  electrolyte 
at  the  higher  temperatures. 

The  Rp  data  for  603XX  are  given  in  Fig.  9.  On  the  basis 
of  their  shapes,  together  with  the  positive  coefficients  in 
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Fig.  7.  Rp  vs.  cycle  count  for  606XX  (XX  =  temperature  on  graph). 
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Fig.  8.  (a)  Rp  vs.  cycle  count  for  803XX  (XX  =  temperature  on  graph);  (b) 
enlargement  of  a  section  of  Fig.  8a  to  show  the  detail. 


Table  3,  these  curves  may  be  indicative  of  a  system  that  is 
mass-transport  limited,  causing  the  concentration  polariza¬ 
tion  to  increase  according  to  some  simple  function  of  cycle 
count. 

The  situation  is  more  complex  at  606XX.  The  curves  in 
Fig.  7  seem  to  asymptotically  reach  a  constant  value.  We 
found  that  the  Rp  data  could  be  expressed  as  an  equilibration 
process  plus  a  degradation  reaction.  This  is  shown  in  Table  3 
by  means  of  two  rows  of  coefficients  for  each  temperature. 
The  upper  row  with  the  negative  coefficient  is  the  equilibra¬ 
tion  process  and  the  bottom,  the  degradation.  The  degrada¬ 
tion  rate  constants  were  calculated  by  subtracting  the  best-fit 
curve  from  the  Rp  data.  Fig.  10  shows  a  typical  plot  of  the  Rp 
data,  the  equilibration  process,  and  the  degradation  reaction. 
It  is  interesting  to  note  that  all  the  equilibration  processes 
follow  t  ^  kinetics,  as  do  most  of  the  degradation  reactions. 
As  indicated  by  the  t  ^  time  dependence,  they  seem  to  be 
diffusion-limited. 

The  remaining  data  from  60%  SOC  (609XX)  and  most 
of  the  80%  SOC  follow  the  same  basic  trend  with  cycle 
count.  The  rate  constants  are  positive,  indicating  that  the 
concentration  polarization  increases  with  time.  The  80960 
data  indicate  that  an  equilibrium  process  may  be  occurring; 
the  rate  constant  is  negative.  However,  we  could  not 
deduce  the  degradation  rate  equation  with  any  degree  of 
confidence. 
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Table  3 


Coefficients  for  Rp  vs.  cycle  count  equations  for  60%  SOC 


SOC  (%) 

Temperature  (°C) 

t 

ln(0 

S'2 

Constant 

r2 

3 

40 

None 

None 

None 

4.18  x  10~la 

50 

7.71  x  1CT3 

3.00  x  10_1 

0.90 

60 

3.37  x  1(T4 

2.92  x  10~‘ 

0.93 

70 

1.76  x  10“6 

2.83  x  10~‘ 

0.91 

6 

40 

-8.31  X  1(T4 

5.23  x  10~2 

0.98 

8.48  x  10“6 

-1.58  x  10~2 

1.00 

50 

-5.27  x  10“4 

4.32  x  10~2 

0.96 

3.86  x  1(T4 

-1.68  x  10“2 

0.95 

60 

-1.16  x  10“4 

2.58  x  10~2 

0.99 

1.17  x  1(T4 

-9.50  x  10~3 

0.99 

70 

-7.29  x  1(T5 

1.89  x  10“2 

0.97 

1.73  x  1(T4 

-9.50  x  10~3 

0.93 

99 

40 

9.17  x  1(T7 

1.89  x  10~2 

0.93 

50 

9.04  x  10-5 

1.57  x  10“2 

0.97 

60 

7.41  x  10“5 

1.34  x  10~2 

0.94 

70 

3.84  x  1CT7 

1.33  x  10“2 

0.96 

a  Considerable  scatter  in  data. 

Table  4 

Coefficients  for  Rp  vs. 

time  equations  for  80%  SOC 

SOC  (%) 

Temperature  (°C) 

ln(r)  r1/2 

Constant 

r2 

3 

40 

6.50  x  10“4 

3.65  x 

10“2 

0.95 

50a 

7.25  x  10"5 

7.00  x 

10~2 

0.91 

60 

1.35  x  10-3 

2.58  x 

10_1 

0.97 

70 

1.03  x  10-4 

4.98  x 

10“2 

0.97 

6 

40 

1.13  x  10~3 

-5.37  x 

10_1 

0.99 

50 

5.08  x  10“‘ 

-7.84  x 

10~2 

1.00 

60 

1.76  x  10~4 

2.08  x 

10“2 

0.94 

70 

8.82  x  10“4 

3.99  x 

10_1 

0.98 

9 

40 

No  trend  obvious  due  to  scatter  in  data 

50 

2.34  x  10~4 

4.73  x 

10“3 

0.99 

60 

-3.19  x  10~3 

9.08  x 

10-2 

0.93 

70b 

4.31  x  10~5 

8.09  x 

10~3 

0.96 

aThis  could  also  be  represented  as  Rp  =  1.80  x  10  3  ln(f)  +  5.98  x  10  2,  r2  =  0.90. 
b  This  could  also  be  represented  as  Rp  =  8.38  x  10-4  In (t)  +  3.74  x  10-3,  r2  =  0.97. 
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Fig.  9.  Rp  vs.  cycle  count  for  603XX  (XX  =  temperature),  where  the  60340  data  have  been  omitted. 
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Cycle  Count 

Fig.  10.  Plot  of  Rp  vs.  cycle  count  for  the  data,  an  equilibration  process 
and  a  degradation  reaction. 

4.  Conclusions 

From  the  voltage  response  of  the  cells  to  the  life-cycle 
profile  at  40,  50,  60,  or  70  °C,  we  have  separated  the  overall 
impedance  rise  into  two  simpler  terms,  R0  and  Rp.  The  R0 
data  tend  to  follow  expected  temperature-dependent  trends 
(40  >50>60>70°C).  The  Rp  data  trends  show  that  Rp 
can  either  decrease  or  increase  asymptotically  with  time,  but 
the  overall  temperature-dependent  behavior  is  similar  to  that 
of  R0. 

Specifically,  at  60%  SOC,  R0  can  be  expressed  as  an 
equation  of  the  form,  aln(t)  +  bt +  c.  This  form  of  the 
equation  implies  that  there  is  an  initial  rapid  reaction  (e.g. 

SEI  formation)  that  slows  down  with  time  into  diffusion- 

1 10 

limited  kinetics.  The  coefficients  of  both  the  ln(f)  and  t 
terms  from  the  3  and  6%  ASOC  data  display  Arrhenius-like 
behavior.  At  80%  SOC,  the  equation  for  R0  is  simpler, 
at +  b ,  with  the  coefficients  of  t1  2  displaying  Arrhe¬ 
nius-like  behavior  in  the  temperature  range  of  40  to  60  °C. 

For  Rp  at  60%  SOC,  two  processes  may  be  occurring: 
an  equilibration  and  a  degradation.  The  former  follows 
at  +  b  or  at1!2  +  b  kinetics,  while  the  latter,  follows 
at1'2  +  b  only.  Only  the  data  from  6%  ASOC  displayed 
clear  evidence  of  the  degradation  process.  At  80%  SOC,  the 
equilibration  process  followed  either  a  ln(f)  -{-box  at1!2  +  b 
kinetics.  There  was  no  clear  evidence  of  degradation  pro¬ 
cesses. 

We  have  illustrated  the  types  of  processes  occurring  in 
one  Li-ion  cell  chemistry.  Based  on  the  initial  rates,  the 


processes  are  complex  and  inter.  The  R0  term  dominates  the 
dependent  observable  cell  impedance,  but  Rp  adds  a  non¬ 
trivial  contribution.  Together,  they  describe  how  the  cell 
performs  and  how  its  performance  changes  with  time. 
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